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We investigate vortex pinning in thin superconducting films with a square array of rectangular 
submicron holes {"antidots^'). Two types of antidots are considered: antidots fully perforating the 
superconducting film, and "blind antidots", holes that perforate the film only up to a certain depth. 
In both systems, we observe a distinct anisotropy in the pinning properties, reflected in the critical 
current /c, depending on the direction of the applied electrical current: parallel to the long side 
of the antidots or perpendicular to it. Although the mechanism responsible for the effect is very 
different in the two systems, they both show a higher critical current and a sharper IV-transition 
when the current is applied along the long side of the rectangular antidots. 



I. INTRODUCTION 

Type-II superconductors (SC's) with nano-engineered 
artificial pinning arrays are good candidates to study the 
fundamentals of vortex pinning, since, within the limits 
of the lithographic process used for their fabrication, the 
pinning centers can be shaped and positioned in the SC at 
will. Often used systems in that respect are, for example, 
periodic arrays of submicron holes (antidots) ^'^'^'^'^ or 
magnetic dots, placed underneath or on top of the SC 

In SC's with a periodic pinning array (PPA), so-called 
matching effects^ occur, at specific magnetic fields gen- 
erating a number of vortices which "matches" the num- 
ber of available pinning sites. At these integer Hn and 
fractional Hp/q matching fields, the vortices form reg- 
ular geometrical patterns, commensurate with the pin- 
ning array. This strongly reduces the vortex mobility 
and consequently increases the critical current /c. These 
commensurability effects have been recently intensively 
studied for square or triangular arrays of antidotsSiii or 
magnetic dots^^"^ 



In case of a square PPA, the equivalence between 
the two important in-plane directions, [10] and [01] {x- 
and 7/-axis), is conserved if square or cylindrical anti- 
dots (or out-of-plane magnetic dots) are used. However, 
as shown by recent numerical simulations, anisotropic 
pinning properties can be introduced in an otherwise 
isotropic superconducting film by using for example a 
rectangular array of isotropic pinning sites"*^"*^. In the 
present work we study, by means of experiments and nu- 
merical simulations, the anisotropy in the pinning prop- 
erties of a SC film caused by a square array of rectangular 
antidots and blind antidots. 

This paper is organized as follows. First, we con- 
sider rectangular antidots that perforate the supercon- 
ductor completely. For this system, we present exper- 



imental results and use numerical simulations to gain 
more insight in the mechanisms responsible for the ob- 
served anisotropy in the vortex mobility and dynamics. 
In a second part, we investigate rectangular antidots that 
do not perforate the superconducting film completely 
{^^ blind antidots^^). Here, we will present only numeri- 
cal simulations. Experiments on a square array of blind 
rectangular antidots have, to our knowledge, not been 
performed sofar. 



II. SQUARE ARRAY OF RECTANGULAR 
ANTIDOTS 



A. Electrical transport measurements 

1. Experimental details 

We patterned a SC Pb film in a 5 x 5 mm^ cross- 
shaped geometry (see Fig. QJa)) to allow electrical trans- 
port measurements in two perpendicular current direc- 
tions. The central part of the cross consists of two 
300 fim wide strips containing the square array of rect- 
angular antidots (see dark gray area in Fig. ^a)). In 
both strips, the long side of the antidots points in the 
^/-direction. This pattern was prepared by electron- 
beam lithography in a polymethyl metacrylate/methyl 
metacrylate (PMMA/MMA) resist bilayer covering the 
Si02 substrate. A Ge(20 A)/Pb(1500 A)/Ge(300 A) film 
was then electron-beam evaporated onto this mask while 
keeping the substrate at 77 K. After a liftoff process in 
warm acetone, the structure was covered with a thick 
insulating Ge(1000 A) layer for protection against oxida- 
tion. 

Figure ^c) shows an atomic force microscopy (AFM) 
topograph of a 6 x 6 jam^ area of the square antidot lattice 
with a period of 1.5 /im. We see that the antidots have a 
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FIG. 1: Layout of the Pb film with a square array of rectan- 
gular antidots. (a) Cross-shaped geometry of the sample to 
allow for transport measurements in the x- and ^-direction, 
(b) Schematic representation of a unit cell (1.5 x 1.5 /xm^) of 
the antidot array, (c) Atomic force micrograph of a 6 x 6 fim^ 
area of the antidot lattice. 

rectangular shape (0.6 x 1.1 /im^) with rounded corners. 
As shown in the schematic representation of a unit cell of 
the array in Fig.mb), the superconducting paths between 
the antidots are 0.4 fim and 0.9 /im wide, for the x- 
and ^/-direction, respectively. The rms roughness value 
of the Ge/Pb/Ge film, in between the antidots, is less 
than 15 A on a 1 /im^ area. 

The transport measurements were performed in a ^He 
cryostat equipped with a 9 T superconducting magnet. 
The magnetic field was applied perpendicular to the film 
surface. Four-point resistivity measurements were car- 
ried out using an ac bridge at a frequency of 19 Hz and a 
current of 10 jaA. We obtained the superconducting criti- 
cal temperature (Tc = 7.26 K) using a resistance criterion 
of 10 % of the normal state resistance Rn- The ratio of 
the Rn values for the x- and ^/-direction is 

R^/Rl = 1.06 Q/0A6 ft = 2.3, 

which is in excellent agreement with the value 2.25 that 
can be expected from geometrical considerations. To 
determine the superconducting coherence length ^(0), 
we measured the linear Tc{H) phase boundary of a co- 
evaporated reference film without any in-plane nano- 
structuring. From the Tc{H) slope, and^^ 

M0^fc2 - ^^^^ - \^~tJ' 

we find ^(0) = 51 nm. Using the dirty limit (^ < ^o) ex- 
pression ^(0) = 0.865^^ and the BCS coherence length 



for Pb, ^0= 83 nmi^, we determined the elastic mean free 
path ^ = 42 nm. We can derive the penetration depth 
A(0) = 34 nm by means of the dirty limit expression 
A(0) = 0.66Alv^o/^, using Al=37 nm as the London 
penetration depth^^. Taking into account that the pres- 
ence of antidots in a superconducting film has the ten- 
dency to increase the penetration depth in the following 
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wa}«s: 

A(0) = /^^^ = 53 nm , (2) 

where Sa is the area occupied by the antidots and St is 
the total area of the film, we obtain a Ginzburg-Landau 
(GL) parameter n = A(0)/^(0) ^ 1 > 1/^2. We there- 
fore conclude that the patterned film with the array of 
rectangular antidots is a type-II superconductor. 



2. Results 

In Fig. m we show the critical current versus field 
curves Ic{H), normalized to its value at zero field, Ico = 
Ic{H = 0), at two temperatures {T/Tc = 0.995 in (a) 
and 0.992 in (b)). We have used a voltage criterion of 
Vcrit = 100 Rn /jM/ft^ in order to be able to make a 
comparison between the two current directions with dif- 
ferent normal state resistances R^ and i?^. The criti- 
cal current density at zero field measured along the x- 
direction is I^^ = 4. A ■ 10^^ at T/T^ = 0.995 and 
I^^ = 9.7 • 10^^ at T/Tc = 0.992. For a current along 
the TV-direction, these values are IL = 1.3 • lO^A at 
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T/Tc = 0.995 and = 2.0 • 10^^ at T/Tc = 0.992. 
The field axis is given in units of the first matching field 
i^i, the field at which the density of ^o-vortices equals 
the density of antidots: fioHi = ^^/d? = 9.2 Oe, with 
d = 1.5 /im the period of the antidot lattice and 0o = ^ 
the superconducting flux quantum. 

Due to the presence of the antidot lattice, both the 
Icx{H) and Icy{H) data show pronounced maxima at 
the integer matching fields Hi and H2 up to the same 
critical current value. The differences between the two 
current directions appear in the field ranges between the 
integer matching fields. For those field intervals, the crit- 
ical current Icy{H) (filled symbols in Fig. [2|) is consider- 
ably enhanced compared to Icx{H) (open symbols). This 
increase of Icy is accompanied by the complete suppres- 
sion of the rational matching peaks at i^p/g {p and q 
integers). For Icx{H)^ on the other hand, the rational 
matching peaks at ^1/3, ^2/3 (panel (a)), and also 

at 3/2 Hi (panel (b)) are clearly revealed. 

To investigate the origin of this qualitative difference in 
the behavior of the critical currents along the two direc- 
tions, we have a closer look at the V{Ix) (open symbols) 
and V{Iy) (solid lines) curves for some selected field val- 
ues at T/Tc = 0.995 (Fig. To make a comparison 
possible, the voltage axis is rescaled to Rnlco^ and the 
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FIG. 2: Normalized critical currents lex and Icy as a function 
of normalized magnetic field H/Hi for the sample shown in 
Fig. measured with a current in the x- (open symbols) and 
the y-direction (filled symbols), respectively. The Ic{H) data 
are presented for two temperatures: (a) T/Tc = 0.995 and (b) 
T/Tc = 0.992. 



current axis to Ico- The dotted horizontal line indicates 
the voltage criterium that was used to define the critical 
current shown in Fig. [21 

At = and H = Hi, the V{I) transitions are very 
sharp and independent of the current direction. The 
V{Ix) and V{Iy) curves in Fig. [31 are therefore almost 
indistinguishable. At H = Hi, the appearance of a low- 
voltage tail in the V{I) curve can be observed due to a 
small deviation from the first matching field. These V{I) 
curves have a shape that is typical for a regular pinning 
arra y -^^1-^^ as it was e.g. also observed in SC films with a 
periodically modulated thickness^^. 

For H = 0.4 Hi, the current at which the film com- 
pletely reaches the normal state is the same for the two 
current directions, but a large tail in the V{Ix) curve 
(open symbols) is present, indicating the dissipative mo- 
tion of vortices in a direction perpendicular to the cur- 



rent. Consequently, the V{Ix) transition is very broad 
(transition width AI = 0.8 Ico) compared to a much 
smaller broadening for V{Iy) (A/ = 0.2 Ico)- This be- 
havior is typical for every magnetic field in between the 
matching fields. 

The field dependence of the transition width AI can be 
more adequately shown by presenting the V{Ix-,H) and 
V{Iy, H) curves in a contour plot (Fig.^IJa) and (b)). We 
see that, at = and H = Hi, the V{I) transitions 
are equally sharp for both current directions. In between 
the integer matching fields, a substantial broadening of 
the V{Ix) transition takes place (Fig. El^a)). At the ratio- 
nal matching fields i^i/s, ^^1/2, and i?2/3: the transitions 
regain part of their sharpness. For / || y (Fig. ^b)), all 
V{Iy) curves, both at and in between the integer match- 
ing fields, have a sharp transition. 

Summarizing the experimental observations, we have 
found a strong anisotropy in the critical current Ic{H) 
and the V{I) characteristics of a film with a square array 
of rectangular antidots. At the integer matching fields 
Hn, both the critical current Ic{Hn) and the V{I,Hn) 
curve do not depend on the direction of the applied cur- 
rent - / II X or / II y. However, in between the integer 
matching fields, we find broad V{Ix) transitions, a low 
lex J and rational matching features in lex {H). For / II y, 
on the other hand, we see sharp V{Iy) transitions at ev- 
ery magnetic field, an overall high Icy, and no sign of 
rational matching features in Icy{H). 




0,0 05 LO 



I/I 

CO 

FIG. 3: Normalized V{I) curves at selected magnetic field 
values at the same temperature (T/Tc — 0.995) as the mea- 
surements in Fig. |2Ia). The full lines are the data for / || y, 
the open symbols for / || x. The horizontal dotted line indi- 
cates the voltage criterium used to define the critical current. 
For if = and H = Hi, the curves for the two current 
directions overlap almost completely. 
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FIG. 4: Contour plots (a) V(Ia:,H) and (b) V{Iy, H) at T/Tc 
= 0.995. The gray scale indicates the normalized voltage 
Vn — p • The evolution of the Vil) transition width as a 
function of magnetic field can easily be seen from this plot. 



S. Discussion 

The differences between Icx{H) and Icy{H) are due to 
an anisotropy in the vortex mobility in the sample along 
the two in-plane directions. Neglecting the thermal noise 
force, the velocity v of a vortex in a PPA is determined by 
the superposition of three forces: the vortex- vortex inter- 
action Fyy, the vortex- ant idot pinning force Fp and the 
Lorentz force Fi,, directed perpendicular to the applied 
current. The critical current Ic of the system is propor- 
tional to the Lorentz force F^ which is needed to induce 
a threshold average vortex velocity, and therefore a cer- 
tain critical voltage. To reach this voltage, a noticeable 
fraction of the vortices must be depinned by the joint ef- 
fect of the Lorentz ¥l and the vortex- vortex interaction 
Fvv forces. The balance between these forces, and con- 
sequently depends strongly on the applied magnetic 
field, since the latter determines the density of vortices 
in the sample. 



When, for example, all pinning sites in a square PPA 
are occupied by a vortex (first matching field i^i), the 
vortex lattice has such a high symmetry, that all vortex- 
vortex interactions between the vortices trapped in the 
pinning sites cancel out {Fyv = 0). Therefore, at this 
field, the Lorentz force Fl which is needed for depinning, 
is only determined by the pinning force. This implies that 
the critical current at the first matching field /c(i^i) is a 
measure of the single site pinning force of the individual 
antidots. Once the Lorentz force exceeds the pinning 
force, all vortices leave their pinning site simultaneously, 
resulting in a sharp V{I) transitions^. 

Since the Icx{Hi) and Icy {Hi) values and the 
V{Ix,Hi) and V{Iy,Hi) curves (Figs. [3 and coincide, 
we conclude that the pinning force exerted by a rectan- 
gular antidot on a single pinned ^o-vortex is isotropic 
along the two symmetry- axes of the rectangular antidot. 
This experimental observation is in agreement with pre- 
dictions by Buzdin and Daumens^^ for an elliptic antidot, 
where the pinning force is only expected to be anisotropic 
when the ellipse is extremely elongated. 

The observed anisotropy in pinning properties can 
therefore not be attributed to an anisotropic single site 
pinning force of the rectangular antidots themselves, but 
rather to their arrangement in the array. In other words, 
the observed anisotropy in the pinning properties should 
be associated with an anisotropic vortex-vortex interac- 
tion. Indeed, since the SC strands between the anti- 
dots, where the screening currents of the vortices have 
to fiow, are much thinner between the adjacent antidots 
in the ^/-direction than in the x-direction, we expect the 
^/-component of the vortex- vortex interaction forces to be 
considerably larger than their x-component. 

To understand how the anisotropic vortex-vortex in- 
teraction can give rise to the observed anisotropy in the 
pinning properties in between the matching fields, we 
examine the role of the vortex- vortex interaction in the 
depinning process of vortices in a regular pinning array. 

We have discussed earlier that at a matching field, the 
vortex-vortex interaction vectors are canceled out and 
depinning is governed by the single site pinning force. 
When the applied field is now slightly detuned from a 
matching field (integer or rational), the vortex lattice 
contains defects due to the incommensurability of the 
vortex array and the PPA. The vortex rows, parallel 
to Fl, without such defects are called "commensurate" 
rows. In these commensurate vortex rows, the compo- 
nent of Fyy parallel to Fi, will be zero for all vortices 
in the row. On the other hand, in vortex rows contain- 
ing defects, the vortex- vortex interaction forces will not 
cancel out at all (Fyv 7^ 0). This Fyv component in 
the direction of the Lorentz force Fl will assist to the 
depinning of these "incommensurate" rows, i.e. depin- 
ning occurs when the sum of this Fyy component and 
Fl overcomes the single site pinning force. The com- 
mensurate rows stay pinned up to a higher Lorentz force 
Fl—. If the vortex- vortex interaction is sufficiently large, 
the critical current Ic{H) will therefore be substantially 
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lower right before and after the matching fields, leading 
to pronounced matching peaks (integer and rational). 

When the current is applied in the x-direction 
(Fl II y), the Icx{H) curve shows pronounced maxima 
at the rational matching fields i^i/2, ^3/2? ^1/3 and 
H2/3 and a much lower value in between these rational 
matching fields (see Fig. This is a clear sign of a 
strong vortex-vortex interaction Fyv component along 
the direction of the Lorentz force F^,, in this case the 
^/-direction. 

When a current is applied in the ^/-direction (F^ || x), 
we see a complete disappearance of the rational matching 
peaks and, instead, an overall high critical current Icy{H) 
(Fig. El). This indicates that the Fyy components along 
the direction of the Lorentz force (in this case the x- 
direction) are much weaker than in the case / || x. The 
vortex- vortex interaction is apparently not able to pro- 
vide the long-range order needed to generate the sparse 
geometrical patterns at the rational matching fields i^p/g. 
This implies that no commensurate rows in the direction 
of Fl are formed, and all vortex rows are qualitatively 
equal (all are "incommensurate"). All rows will depin ap- 
proximately at the same driving current. This explains 
the sharp V{Iy) transitions for all magnetic fields (see 
Fig. ^b)). The integer matching peaks are still present 
in Icy{H)^ since the commensurate vortex patterns at the 
integer matching fields Hn can be achieved even with a 
small interaction force present, due to a smaller vortex 
separation. 

Our critical current results are in agreement with what 
was found in phase boundary measurements on square SC 
networks with a different wire width in the two perpen- 
dicular directions'^. In these systems, Tc{H) depends on 
the direction of the current in a similar way as in our 
Ic{H) data. This is not surprising since there are many 
similarities between SC films with an antidot lattice and 
SC wire networks, especially at temperatures close to Tc. 
Using the temperature to tune the coherence length ^(T), 
one can, in fact observe a weakly coupled network behav- 
ior (^(T) ^ strips between the antidots) in any antidot 
array. 

Summarizing this part of the discussion, our main 
experimental observations can be explained by (i) an 
isotropic single site pinning force of the rectangular an- 
tidots, combined with (ii) an anisotropic vortex- vortex 
interaction between the vortices trapped inside the anti- 
dots. 



B. Numerical simulations 

We have used molecular dynamics simulations to con- 
firm that a model system with the features (i) and (ii) 
mentioned above, shows indeed the anisotropy in the crit- 
ical current that was observed in the experiments. 



1. Model 

We consider a two-dimensional system with periodic 
boundary conditions containing a square array (period 
d) of circular (i.e. isotropic) pinning sites, but with an 
anisotropic vortex-vortex interaction. The overdamped 
equation of motion for a vortex i (Eq.Q) is used to cal- 
culate the average vortex velocity as a function of the 
applied Lorentz force F^, and to trace the vortex trajec- 
tories in the pinning array: 



r]Wi = Fl + Fyy(Yi) + Fp(ri), 



(3) 



where 77 is the viscosity coefficient and taken to be unity. 
The driving force acting on the vortices is the Lorentz 
force, = J x (/)o, where J is the applied current. 

The anisotropic repulsive vortex-vortex interaction 
force between two vortices with separation r^j is de- 
scribed by the modified Bessel function ^-,^2i2iSii2&: 



Fyy (r^ 



where r^.- = t^-^^i ^md N denotes the number of vor- 

J I Yi Yj I 

tices in the sample. Fyy^fo expresses the intensity of the 
vortex- vortex interaction force, with /o 



vvO 



(4) 



87r2A3 



. The 



cut-off for the vortex- vortex interaction lies at 6 d in both 
the X- and the ^/-direction. The vortex- vortex interaction 
is made anisotropic by choosing Xy = 2Xx = d. All forces 
are expressed in units of /q. 

The rectangular antidots of the experiment are mod- 
eled as potential wells with an isotropic attractive pin- 
ning force ¥p given by: 



Fp{r) = -Fpofo -—exp 
k=i 



(5) 



Here, r is the distance between the vortex and the k^^ 
pinning site. 

Calculations are performed on a square sample (side 
8 d) with periodic boundary conditions in both x- and 
^-directions. The same initial vortex configuration, ob- 
tained from an annealing course^^, is used for the hys- 
teretic calculation of the {v){Fl) current- voltage charac- 
teristics along the x- and the ^-direction. 

We solve the equation of motion (Eq.Q) by taking 
discrete time steps At, using the following fixed lengths 
and forces: Rp = 0.7, = 6, Fyyo = 0.2, Fpo = 1. The 
average vortex velocity as a function of applied driving 
force {v){Fl) is obtained by increasing Fl in small steps 
AFl = 0.001, starting from Fl = 0. At each F^, we 
neglect the first 2000 time steps, and average the next 
8000 steps to calculate the average velocity of all the 
vortices in the system. 
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2. Results 

The obtained average vortex velocity as a function 
of the apphed Lorentz force for H/Hi = 0.58, i.e. a 
magnetic field in between matching fields, is shown in 
Fig.El We will compare these calculated {v){Fl) curves 
at H/Hi = 0.58 with the experimental V{I) curves taken 
at H/Hi =0.4 (Fig. IS)), which are typical for all magnetic 
fields in between matching fields. The two main features 
observed in the experimental curves are qualitatively re- 
produced. First, the onset of the vortex motion occurs 
at lower driving force for / || x (F/, || y) (open symbols) 
than for / || y {F^ \\ x) (filled symbols). The critical 
depinning force in the ^/-direction lies at Fl ~ 0.315, 
while it is Fl ^ 0.37 when the force is applied in the 
x-direction. Therefore, the critical current for / || x is 
lower than for / || y. Secondly, the V{I) transition is 
clearly much broader for / || x than for / || y, with a 
gradually increasing average vortex velocity. 
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FIG. 5: Average vortex velocity as a function of the Lorentz 
force Fl at H/Hi = 0.58. The open (filled) symbols show the 
result for a current along the x- (y-) axis. Arrows a and c (b 
and d) indicate the points at which the vortex trajectories for 
/ II y (/ II x) were calculated, shown in Fig. El 

The origin of the different {v){Fl) behavior for the two 
directions of the Lorentz force can be found by examin- 
ing the trajectories of the depinned vortices (Fig. E|) at 
two Lorentz force values, indicated by arrows in Fig. [51 
A Lorentz force of 0.32, applied in the x-direction (arrow 
a in Fig. 13) is not sufficient to depin any vortex rows (see 
Fig. ini^a)). When the same Fl = 0.32 is applied in the 
^/-direction (arrow b in Fig.EJ, a partial depinning takes 
place, where the rows with little order in the ^/-direction 
(incommensurate rows) are the first ones to start moving 
(see Fig. jS^b)). The commensurate rows (the first, sec- 
ond, and forth row, counting from the left) remain pinned 
(see also Section FlI A 3|) . At a Lorentz force Fl = 0.38, 
applied along the x-direction (arrow c in Fig. [5|), all the 
vortex rows have started to move (Fig. jS^c)). This cor- 



responds to the sharp increase in (v) that can be found 
in the {v){Fl) curve (filled symbols in Fig. For the 
same Fl = 0.38 applied in the ^/-direction (arrow d in 
Fig. 0), the vortices in the incommensurate rows move 
with a higher velocity (Fig. jS^d)). When an even higher 
Lorentz force is applied (e.g. Fl = 0.44), all vortices 
finally move for both Fl directions (not shown). 



3. Discussion 

The results of these molecular dynamics simulations 
show that an array of isotropic pinning sites^ combined 
with an anisotropic vortex-vortex interaction along the 
two principal axes of the PPA, gives rise to the same 
phenomena as observed in our pinning experiments on 
a Pb film with a square array of rectangular antidots. 
The low critical current and broad V{I) transitions in 
between matching fields for / || x (F^ || y) are found to 
be due to the motion of the incommensurate vortex rows 
along the direction of the strong vortex- vortex interac- 
tion, in this case the ^/-direction. For / || y {Fl \\ x), 
the simultaneous vortex depinning at a relatively high 
driving force is responsible for the high critical current 
and the sharp V{I) transitions. This result strengthens 
our interpretation that in our experiment, the rectangu- 
lar antidots provide an isotropic pinning potential but 
induce an anisotropy in the vortex-vortex interaction. 
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FIG. 6: Vortex trajectories (lines) for H/Hi = 0.58 at Fl — 
0.32 for a current in the y- (a) and the x-direction (b). Open 
circles represent pinning sites, black dots represent vortices. 
Idem at Fl — 0.38 for the current in the y- (c) and x-direction 
(d). The labels (a) to (d) correspond to the arrows in Fig.|51 

Due to the anisotropic nature of the vortex- vortex in- 
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teraction, the initial vortex lattice found in the simu- 
lation by an annealing course shows more order in the 
^/-direction (stronger vortex-vortex interaction) than in 
the x-direction (weaker vortex-vortex interaction) (see 
Fig. Ela)). Although we believe that this is also what 
happens in the experiment, no definitive statement about 
the experimental vortex configuration can be made, how- 
ever, at this stage, before a direct experimental study of 
the vortex patterns in these anisotropic pinning arrays, 
e.g. by scanning Hall probe2Si2i or Lorentz^^ microscopy. 

It is interesting to compare our results with the molec- 
ular dynamics simulations on a rectangular array of 
isotropic pinning sites, recently performed by Reichhardt 
et al}^ . In that case, the anisotropy in Fyy between 
pinned vortices is induced by a different spacing in the x- 
and ^/-direction of the pinning sites. A similar anisotropy 
in Ic{H) as in our measurements and calculations is 
indeed found when no interstitial vortices are present 
{H < Hi). When the Lorentz force is in the di- 
rection of the weak vortex- vortex interaction (in the case 
considered in Refii along the direction with a large lat- 
tice spacing) the absence of rational matching and an 
overall high Ic{H) is found. 



Moreover, a vortex trapped inside a blind antidot can 
move around freely within the boundaries of the pinning 
center. In a film with antidots, on the other hand, the 
supercurrents constituting the vortex are forced to fiow 
at the outer edge of the antidot, confining the vortex in 
space. 



A. Model 

Taking into account the properties of a blind antidot, 
mentioned above, we will use the following model to simu- 
late the vortex motion in a superconductor with an array 
of blind antidots. 

Again, we use a square 2D system (size 8d) with peri- 
odic boundary conditions in the x- and ^/-direction. The 
system contains a square array (period d) of rectangu- 
lar pinning sites (sides Rx = OAd and Ry = 0.77d). 
The vortex velocities are calculated using the normalized 
overdamped equation of motion Eq. The repulsive 
vortex- vortex interaction is described by Eq. using 
an isotropic penetration depth X = d and Fyyo = 0.2. 



C. Summary 

We have measured Ic{H) and V{I,H) curves of a SC 
film with a square array of rectangular antidots for two 
directions of the applied current. We find an overall 
high Ic{H) with integer matching peaks but no ratio- 
nal matching features when the current is applied par- 
allel to the long side of the antidots. This is a clear 
advantage compared to isotropic pinning arrays, where 
a large suppression of Ic{H) is seen in between the ra- 
tional matching fields. We attribute this effect to the 
anisotropic vortex-vortex interaction^ which is stronger 
along the long side of the antidots. Molecular dynamics 
simulations on a square array of isotropic pinning sites, 
combined with an anisotropic vortex-vortex interaction, 
indeed show the same characteristic anisotropic features 
as observed in our measurements. 



III. SQUARE ARRAY OF RECTANGULAR 
BLIND ANTIDOTS: NUMERICAL 
SIMULATIONS 

In this Section, we examine the vortex mobility in a 
superconductor with a square array of rectangular blind 
antidots by means of molecular dynamics simulations. 
Blind antidots are holes that do not perforate the super- 
conductor completely, i.e. the bottom of a blind antidot 
always contains superconducting material. Therefore, if 
more than one vortex is trapped inside a blind antidot, 
the vortices remain individual repulsive entities, without 
merging into one multiple-quantum vortex^'^^, as in the 
case of an antidot. 




FIG. 7: 3D plot of the pinning potential used to model a 
square array of rectangular blind antidots. 

The pinning force of a rectangular pinning center on 
the vortex can be defined as^^: 



Fp(r) 



-i^po/o V -—exp 



k=l 



Rv 



(6) 



when the vortex is outside the k^^ rectangular pinning 
center, and Fp(r) = 0, if the vortex is inside the rectan- 
gular area of the pinning site. This means that no pinning 
force is exerted on the vortex, once it is trapped inside 
the blind antidot. In this expression, r is the distance 
between the vortex and the edge of the k^^ rectangular 
pinning center. We have chosen Rp to be Rp = 0.067d 
and Fpo = 1. In Fig. [3 we show schematically the 3D 
pinning potential that was used for this simulation. 

The current- voltage characteristics {v){Fl) and the 
vortex trajectories for a current along the x- and the 
^/-direction are calculated using the method described in 
Section HTbTI 
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B. Results and discussion 

In Fig. [HI we show the calculated {v){Fl) curves at 
H/Hi = 0.8, for a driving force Fl in the x- (filled sym- 
bols) and the ^/-direction (open symbols). The magnetic 
field is chosen in between matching fields, where no com- 
mensurate vortex configuration can be established. 




FIG. 8: Average vortex velocity as a function of the Lorentz 
force Fl at H/Hi = 0.8. The open (filled) symbols show the 
result for a current along the x- (y-) axis. Arrows a and c {b 
and d) indicate the points at which the vortex trajectories for 
/ II y (/ II x) were calculated, shown in Fig. El 



energetically most favorable vortex configuration looks 
like the configuration at i^3/4, but with some excess vor- 
tices that are placed in the rows with half-filling. Because 
the jumps are assisted by the Lorentz force, we find this 
kind of order in the vortex positions only in the direc- 
tion of the applied Lorentz force (along the x-direction 
in Fig.lHl^a) and along the ^/-direction in Fig.[3fb)). 

A Lorentz force of 0.36 applied in the x-direction (ar- 
row a in Fig. [HI is not sufficient to depin any vortex rows 
(see Fig. IHJa)). When the same Fl = 0.36 is applied in 
the ^/-direction (arrow b in Fig. [H|), a partial depinning 
takes place, where the row with the least order in the y- 
direction (the first vertical row counting from the right in 
Fig.Efb)), is set in motion. The moving vortex row is the 
only one that contains 6 vortices. Accommodating these 
six vortices in a row of eight pinning sites gives rise to a 
lot of strain, making it considerably easier to depin this 
incommensurate row. Most of the other rows contain ei- 
ther eight or four vortices, which can be easily positioned 
in a low energy configuration (commensurate rows). One 
vertical row accommodates five vortices. This row will 
be the next to be depinned when the Lorentz force is in- 
creased (at Fl ~ 0.37), creating the second step in the 
{v){Fl) curve for || y in Fig. [HI (open symbols). 

At a Lorentz force of Fl = 0.4, applied along the x- 
direction, all vortices are set in motion (Fig.Efc)). This 
corresponds to the sharp increase in (v) that is seen in the 
{v){Fl) curve (filled symbols in Fig.IHJ. When a Lorentz 
force of Fl = 0.48 is applied in the ^/-direction (arrow d 



The calculation shows that the onset of vortex motion 
occurs at a lower Lorentz force Fl when F^, || y. For this 
direction of the driving force, the average vortex velocity 
becomes finite at Fi, ^ 0.356. When the Lorentz force is 
applied in the perpendicular direction (F^ || x), vortex 
motion is delayed until Fl ^ 0.386. Consequently, the 
critical current Ic is smaller for / || x than for / || y. 

Moreover, there is a striking difference in the width 
of the depinning transition. When Fl \\ y, we observe 
several steps in the depinning process, leading to a broad 
transition, while for Fi, || x, the {v){Fl) curve shows a 
single sharp jump at Fl ~ 0.386. 

We examine the origin of this qualitatively different 
behavior for the two directions of the Lorentz force by 
plotting the vortex trajectories and positions after 10^ 
time steps (see Fig. 0) at some specific Lorentz force val- 
ues, indicated by arrows in Fig. [HI 

First, we observe a different configuration of the pinned 
vortices in Fig.[3Ja) and (b), although they originate from 
the same annealing course. The difference results from 
vortex jumps that occur after the annealing course, un- 
der the inffuence of the applied Lorentz force. These 
jumps aid the system to obtain its lowest energy configu- 
ration. At a magnetic field of H/Hi = 0.8, it is not pos- 
sible to obtain a completely regular vortex pattern. The 
commensurate configuration, which is closest to 0.8i^i is 
0.75i^i = i^3/4, where one finds consecutively completely 
filled rows and rows with half-filling. At H/Hi = 0.8, the 
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FIG. 9: Vortex trajectories (fines) for H/Hi = 0.8 at Fl = 
0.36 for a current in the y- (a) and x-direction (b). Rectangles 
represent the blind rectangular antidots, black dots represent 
vortices. Idem at Fl = 0.40 for a current in the ^/-direction 
(c), and at Fl — 0.48 for a current in the x-direction. The 
labels (a) to (d) correspond to the arrows in Fig. |H1 
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in Fig.[S|), all vortices participate in the motion. 

In Fig. ini^a), one also finds commensurate horizontal 
rows, containing eight or four vortices, one incommensu- 
rate horizontal row with five and one with six vortices. 
However, these incommensurate rows do not give rise to 
a partial depinning as in the case F^, || y, because the 
strain that is caused by the incommensurability is effec- 
tively released by shifting the vortex positions within the 
rectangular blind antidot. When F^ || x, the long side 
(Ry) of the rectangular blind antidot can be used for this 
purpose. For F^ || y, the vortices are all located at the 
top edge of the rectangular antidots. Consequently, only 
the short side of the antidots {Rx) is available, result- 
ing in large stresses in the incommensurate vertical rows 
which initiate the partial depinning proces. 

We also notice that there is some amount of row switch- 
ing at high driving force for F^ || x due to the smaller 
spacing between the blind antidots in that direction. Al- 
though the plotted trajectory lines in Fig. M,^) do not 
show it directly, we know the switching takes place by 
comparing the number of vortices in a particular hor- 
izontal row at low and at high driving force. On the 
other hand, no row switching is seen when F^ || y. 

We have also examined the effect of the side ratio of 
the rectangular blind antidots Ry/Rx on the vortex de- 
pinning. In Fig. 03 we present {v){Fl) curves (shifted 
vertically for clarity) at H/Hi = 0.8 for various Ry j Rx •, 
smaller than the ratio used for Fig. [HI 

When Fl || y (Fig. EHIb)), the (v)(Fl) curves are al- 
most independent on the ratio Ry/Rx- The transition is 
always a two-step process, involving partial depinning of 
the incommensurate vertical rows, as described above. 

When F/, II X (Fig. QHla)), we find a sharp depinning 
transition, resulting from simultaneous depinning of all 
vortices and a high critical current, if Ry is sufficiently 
large, Ry > 1.8 (see also Fig. |H1 where Ry/Rx = 1.92). 
For Ry/Rx < 1.5, two step transitions in the {v){Fl) 
curves are always seen. The onset of vortex motion also 
appears to be strongly dependent on the ratio Ry/Rx- 
As long as Ry > Rx^ the critical current is higher when 
¥l II X. When Ry < Rx^ the highest critical current is 
found for Fl || y. 

When Fl \\ y, the important length scale determin- 
ing the confinement of the trapped vortices is the width 
Rx of the rectangular blind antidots. Since Rx is kept 
constant in the calculations, we indeed expect to find no 
dependence of the {v){Fl) transitions on Ry/Rx. 

When ¥l \\ x, on the other hand, the initial critical 
depinning force is very sensitive to Ry^ because now this 
parameter determines the confinement of the vortices. 
When Ry/Rx = 0, i.e. Ry = 0, the pinning centers are in 
fact points in the x-direction, which results in a very low 
threshold value for the initial vortex depinning transition. 
When Ry/Rx increases, the vortices trapped in the rect- 
angular blind antidots obtain more freedom to shift their 
positions along the ?/-axis. This gradually improves their 
ability to reduce the strain arising from the vortex- vortex 
interaction in an incommensurate row. In this way, the 
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FIG. 10: Average vortex velocity as a function of the Lorentz 
force Fl at H/Hi = 0.8 with Fl || x (a) and Fl || y (b). 
Rx = OAd is kept fixed and the ratio Ry/Rx is varied from 
0.0 to 1.8, as indicated. 



initial vortex depinning of the incommensurate rows can 
be delayed to higher F^. 



C. Summary 

We have investigated numerically vortex pinning and 
dynamics in a square array of rectangular blind antidots 
for different orientations of the applied current or Lorentz 
force. We demonstrate that the onset of vortex motion, 
and therefore the critical current, depends on the orien- 
tation of the Lorentz force. The critical current is found 
to be higher when the current is applied parallel to the 
long side of the antidots. By shifting their position in- 
side the blind antidot, the vortices are able to delay the 
depinning of the incommensurate vortex rows to higher 
Fl , by reducing the strain arising from the vortex- vortex 
interaction in the incommensurate rows. We have also 
found a strong dependence of the {v){Fl) transition on 
the aspect ratio of the rectangular blind antidots. 
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IV. COMPARISON OF THE ANTIDOT AND 
THE BLIND ANTIDOT SYSTEM 

This paper presents results on the anisotropy in vortex 
pinning for two very different systems: a superconducting 
film with a square array of rectangular antidots, and with 
a square array of rectangular blind antidots. We have 
shown that in both systems the critical current is higher 
and the IV-transition is sharper when the current is sent 
parallel to the long side of the antidots. On the other 
hand, a current along the short side of the rectangular 
antidots yields a lower critical current and a broad IV- 
transition in both systems. 

Although the behavior of the antidot and blind antidot 
system appears to be similar, the mechanisms determin- 
ing this behavior are very different in the two systems. 
In the case of the superconducting film with rectangular 
antidots, the experiments show that the pinning force of 
the rectangular antidots is identical for the two perpen- 
dicular directions. Here, the vortex-vortex interaction is 
found to be anisotropic and responsible for the observed 
anisotropy in the critical current measurements. In the 
case of rectangular blind antidots, the vortex-vortex in- 
teraction is isotropic. Here, the anisotropy in the shape of 
the antidots and the fact that a vortex can move around 
freely within the blind antidot cause the anisotropy in 
the critical current. 



V. CONCLUSION 

We have investigated vortex pinning and dynamics in 
superconducting films with a square array of rectangular 



antidots and blind antidots. As shown by experiments 
(on the antidots system) and numerical calculations (on 
both systems), the critical current is higher when ap- 
plied parallel to the long side of the antidots. For this 
current direction, we also find sharper IV transitions in 
both systems. By examining the vortex trajectories after 
depinning by means of numerical simulations, we were 
able to demonstrate that the origin of the anisotropy is 
very different in the antidot and the blind antidot system. 
Nevertheless, both systems illustrate how patterning of 
an otherwise isotropic superconducting film, can induce 
an important anisotropy in the critical current and IV 
characteristics. 
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